Search for top in lepton + jets in D0 using a topological tag 
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Abstract 

We have searched for production of tt pairs in pp interactions at 1.8 TeV center-of-mass energy at 
the FNAL Tevatron collider. The search assumes standard model decay for top quark into W + b 
quark. We observe in e + jets and [i + jets final states a small, not statistically significant, excess 
above the background estimated by two different methods. The results presented are preliminary. 
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ijj Assuming the Standard Model decay for top quarks into 
a W and a b quark, events with tt pairs will have one 

^ isolated high pt e or /i accompanied by multiple jets 30% 
•i-H of the time. These channels yield considerably more 
^n] events than those with 2 isolated leptons in the final 

^ state, presented in the preceding contribution to these 
" proceedings by S. Wimpenny Q , but suffer from larger 
backgrounds from the channel W+ multi-jets. In this 
paper we discuss two different methods used to estimate 
whether a significant tt production can be observed in 
the D0 data over the IU+jets background. A previous 
analysis of these data, searching for top in the mass 
region 90 to 140 GeV in single and di-lepton final states, 
has been published [Q] and gave a top mass limit of 
131 GeV (95% CL). The present analysis is optimized 
for masses above 130 GeV and should be considered 
preliminary. The two methods used yield a small, but 
not statistically significant, excess. 



2. Event selection 

The D0 detector is described in detail in || and in these 
proceedings by so it will not be discussed here. The 
events for lepton + jets analysis are obtained with 3 
triggers: 



One electromagnetic shower with transverse energy 
(E T ) > 15 GeV, one hadronic jet with E T > 10 GeV 
and calorimeter missing Et (^t 11 ) > 10 GeV 
One electromagnetic shower with Et > 20 GeV and 



if 1 > 20 GeV 



• One fi with px > 8 GeV and one hadronic jet with 
E T > 15 GeV 

Hadronic jets are defined as energy in a cone of R = 
y/ Acj) 2 + A77 2 = 0.5 where 4> is the azimuthal angle in 
radians and rj is the pseudo-rapidity. 

After reconstruction the events are selected further 
by making stricter requirements. For the e+jets sample 
selection requirements on the electron are based on 
shower shape, cluster isolation, shower match with a 
reconstructed track and track ionization cuts to remove 
conversions. Electrons are required to have E T > 20 
GeV and |r; e | < 2.0. Additional requirements on the 
e + jets sample are $J? 1 > 25 GeV, E 3 T et > 15 GeV and 



\r]jet\ < 2.0. Figure |l|a (|l|b) shows the I^t distribution 
for the selected events which pass (fail) the electron cuts. 
We will refer to the events that satisfied the triggers but 
failed the offline electron cuts and have I$ T < 25 GeV 
as the QCD sample. This sample is used to estimate the 
non-W background in the e+jets sample and to test the 
validity of jet multiplicity scaling in section ||. 

Selection requirements on the /i are made to remove 
cosmic rays (by drift chamber timing and no other 
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e + jets 


non-W bckg. 


fi + jets 


non-W bckg. 


> 1 
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94. ± 6. 


303 


48. ± 8. 


> 2 


243 


26 ± 3. 


92 


16. ±3. 


> 3 


35 


5.7 ±0.9 


22 


4.8 ±0.8 


> 4 
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1.1 ±0.3 
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1.3 ±0.3 


> 5 
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0.3 ±0.1 
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0.3 ±0.1 



Table 1. Selected e and fj, + jets events 
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backg. 


W + jets 


tt 


> 1 


1686 


142. ±20.0 


1656 ± 96 


5.7 ± 10.0 


> 2 


335 


50.6 ±7.0 


278.8 ± 17 


5.6 ±9.9 


> 3 


57 


10.3 ± 1.5 


42.2 ± 7.1 


4.5 ±8.2 


> 4 


12 


2.4 ±0.4 


6.6 ±3.8 


3.0 ±4.6 



Table 2. Estimated numbers of W + jets and tt events 



Figure 1. (a) $|? ! distribution for events satisfying all electron 
selection criteria, (b) Tf,?^ 1 distribution for events satisfying the 
electron triggers but failing offline selection criteria. 



collinear /x-track), isolation cuts to remove /i's from 
decays other than W's ( minimum separation between 
fj, and jets with E T > 8 GeV of AR > 0.5 and 
core cone isolation energy < 5 GeV in annular cone 
0.2 < AR < 0.4 around /i track) plus calorimeter dE/dx 
confirmation to reduce accidentals. The fi is required to 
have p T > 15 GeV and < 1.7. In addition to the /z 
the /z+jets sample is required to have at least 1 jet with 
E T > 15 GeV and \r) jet \ < 15 GeV, # T > 20 GeV, and 
Ifi T al > 20 GeV. The trigger and selection requirements 
on the /i+jets sample produces a bias for < 2 jets which 
needs to be corrected when added to e + jets for the jet 
multiplicity analysis in section ||. 

To avoid overlap with the analysis of events in 
which jets are tagged as containing heavy quarks by 
the observation of a non-isolated /x, discussed in the 
contribution to these proceedings by R. Raja /i- 
tagged events have been excluded from the sample and 
their loss accounted for in the acceptance calculation. 

3. Jet Multiplicity Analysis 

The number of events as function of jet multiplicity 
and the estimated non-W background to e + jets and 
fi + jets channels are given in table [j]. The non-W 
background for e + jets is mostly from misidentified 
electrons. In /z+jets events the background comes from 
two sources: Z + jets where one of the fi is not detected 
and QCD processes where the fi from a decay other than 
W satisfies the isolation criteria, the two backgrounds 
are roughly comparable. 

As shown on Fig. 0(a) the e + jets data after 



background subtraction are fit well by a scaling law (as 
expected by 



W + n jets = 



[W+(n- 1) jets) ■ [W + (n - 1) jets)} 
[W+(n- 2) jets)] 



This law seems to be satisfied for any Ei^ threshold 
used for counting jets. The predictions from the 
VECBOS H Montecarlo agree well with the data after 
ISAJET § is used for hadronization and GEANT § for 
simulating the D0 detector. As tt events in channels 
with one isolated lepton are expected to have high 
multiplicities (2 jets from b's, 2 jets from hadronic W 
decay plus additional jets from initial and final state 
radiation) any noticeable tt production should appear 
as an excess of events at high multiplicity, in particular 
for 4 or more jets and for high Et threshold. As shown 
in Fig. ^| the 4th jet (ordered in Et) from high mass top 
is expected to be considerably more energetic than for 
W + jets but no significant violation of the scaling law is 
observed in Fig. |3|(a) . To find how much tt production 
can be accomodated and still satisfy the scaling law we 
can solve using the number of events with 2,3 and 4 jets 
for the number of W and top events: 

V 4 - a 4 • N t = {N 3 - a 3 ■ N t ) 2 /(N 2 - a 2 ■ N t ) 

where TVj is observed number of lepton + jets events with 
>i jets, Nt the total number of top events and the 
corresponding fraction at each multiplicity. The results 
are given in table [|, of particular interest is the total 
number of top events Nt = 5.7 ± 10.0. The numbers for 
> 1 jets are predicted from the numbers obtained using 
the events with 2,3 and 4 jets. 

To check the validity of the scaling law in a 
case with no tt events, we fitted the QCD sample 
to a simple exponential. In Fig. §(b) we show 



Figure 2. Et distributions of 3 rd and 4 th jet from W+ > 3 jets 
Montecarlo (solid) and tt (dotted) with top mass 160 GeV, 
normalized to 500 p6 _1 . 

the jet multiplicities for that sample (after correcting 
for the higher probability as a function of increasing 
multiplicity for having a jet fake an electron). The 
QCD data is clearly well fitted. We also fitted the 
predictions of the VECBOS Montecarlo. All the fits give 
similar values for the ratio of [W + 3 jets]/[W + 2 jets]: 
.151 ± .020 for the data, .160 ± .013 for VECBOS 
Montecarlo, .187 ± .004 for the QCD sample. The fit 
to the VECBOS Montecarlo tends to underestimate 
the predicted number of 4 or more jet events by 15% 
while the QCD fit overestimates by the number of 4 
jet events by 12%. To account for these deviations we 
estimate a 20% systematic error in the validity of the 
scaling law. The fraction (at) of top events expected 
at each multiplicity is obtained from Montecarlo; to 
estimate the uncertainty in this procedure we compared 
2 different Montecarlos (ISAJET and HERWIG @). 
The differences between them were at the level of 
10%. There is an additional 15% uncertainty in the 
top acceptance from the uncertainty in the jet energy 
scale. From the number of top events and the calculated 
acceptance times branching ratio we can estimate a cross 
section for tt production: 

<7u = 6.4 ± 9.8 (stat.) ± 4.0 (sys.) pb 

the sytematic error includes a 12% uncertainty in the 
luminosity. 

This cross section is for masses between 160 and 180 
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Figure 3. Number of events as function of jet multiplicity (a) 
e + jets data after non-W background subtraction and VECBOS 
predictions, boxes indicate uncertainty in prediction, (b) QCD 
sample after correcting for increased probability for triggering 
with increasing multiplicity. 

GeV. As the acceptance falls moderately below 160 GeV 
it should be increased by 25% around 140 GeV. 

4. A, Ht Analysis 

Another method to estimate the number of tt and 
W + jets events in the lepton + jets data is to use event 
shape information to discriminate between them. As the 
ratio of signal/background is small for jet multiplicities 
below 4, this analysis concentrates on events with 4 
or more jets. Two parameters found to be useful 
discriminants are Ht, defined as the scalar Et sum over 
jets with \r]\ < 2.0, and the aplanarity (A), defined as 
3/2 times the smallest eigenvalue of the laboratory 3- 
momentum tensor of reconstructed objects, normalized 
to unit trace. Fig. [| shows the distribution of events 
in this two variables for the QCD sample, VECBOS 
Montecarlo, tt events with mass 180 GeV and the data. 
If we divide the A, Ht into 4 quadrants using axes 



100 200 300 4C0 100 200 300 400 10 15 20 25 30 35 40 45 50 




Figure 4. A vs Ht distributions for events with 4 or more jets: 
(a) QCD sample, (b) VECBOS Montecarlo, (c) tt events, top 
mass 180 GeV, (d) data, the quadrants used in the fit are 
numbered, see text. 



Figure 5. Et distributions for 3 rd and 4 th jet in events with 
> 4 jets starting from VECBOS with 3 jets (dotted lines) and 
from VECBOS with 4 jets (solid lines). 



A = 0.05 and H T = 140 GeV we see that the QCD 
sample and the W + jets Montecarlo sample populate 
those quadrants more or less equally while the tt events 
are concentrated in the quadrant with A > .05 and 
Ht > 140 GeV. In table || we give the expected fractions 
of W + jets , QCD and tt events in each quadrant. To 
estimate the uncertainties in the + jets estimation we 
generated W + jets events in two different ways: one 
was to use VECBOS to generate W + 3 jets events 
(at the parton level) and use ISAJET with those events 
as input to generate events with 4 or more jets. The 
other was to use VECBOS to generate W + 4 jets 
and then proceed as in the previous case. One can see 
from Fig. || that both sets generate quite similar Et 
distributions for the 3rd and 4th jet in events with > 4 
jets. As shown in the table [| the fractions estimated 
with the two different Montecarlo samples differ by less 
than 20%. The fractions for tt were estimated using 
ISAJET, comparison of ISAJET with HERWIG showed 
the differences to be less than 10%. Figure ^ illustrates 
this for tt events generated with ISAJET and HERWIG 
at a top mass of 140 GeV. The number of tt events in 
a given quadrant i {Nli) is given by 

where N is the total number of observed events with 4 or 
more jets, is the fraction of tt expected in quadrant 




Figure 6. A and Ht distributions for tt events with 4 or more 
jets for top mass 140 GeV. Events generated by ISAJET are in 
solid lines, by HERWIG in dotted lines. 



i and / t j is the fraction of N that are tt events. The 
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A > .05 
H T > 140 


e 2 

A > .05 
H T < 140 


A < .05 
Ht < 140 


A < .05 
H T > 140 


VECBOS(4) 
VECBOS(3) 
QCD 


0.21 ± .025 
0.19 ± .04 
0.19 ± .04 


0.27 ± .03 
0.28 ±.05 
0.25 ±.05 


0.21 ± .025 
0.25 ± .05 
0.28 ± .05 


0.31 ± .03 
0.28 ± .05 
0.28 ± .05 


tt 180 


0.60 ± .05 


0.02 ±.01 


0.02 ± .01 


0.36 ± .04 



Table 3. Fractions in > 4 jets 



number of expected background events is then: 



N, 



bkg 



(1 - ftt) ■ N 



where N£ k is the number of background events in 
quadrant i and e\ kg is the expected fraction. Given N 
and the e's one can fit for ] t t using Poisson statistics. 
The results of the fit give f t i — 0.32 ± 0.30 events in 
the sample of data with 4 or more jets. The numbers 
of observed events (in the 4 quadrants of Fig. 4) are 
4, 1, 3, 4 while the predicted numbers are 4.1, 2.4, 
2.0, 3.5. There is a large systematic error (35%) from 
the choice of partitioning the A, Ht plane, this error 
is estimated by moving the axes until one event falls 
into a different quadrant. An additional 20% systematic 
error is estimated for the uncertainty in calculating the 
fractions for each quadrant. The number of tt events 
estimated in the sample with 4 or more jets is thus 
Ntt — 3.8 ± 3.6(stat.) ± 1.5(sys.) which leads to a tt 
production cross section of 

a t t = 8.1 ± 6.9(stat.) ± 3.8(sys.) pb 

additional systematic errors in this estimate are the 
luminosity (12%), and the tt acceptance (20%). This 
value is in good agreement with that obtained with a 
totally different method in section |j. 

The two methods used to estimate a t i can also be 
used as an estimate of the number of background events 
after applying cuts of A > .05 and Ht > 140 GeV which 
improve the signal/background ratio by a factor of 3. 
For method 1 (section ||) we add the fraction of W+ > 4 
jets (predicted from the fit) and the fraction non-W 
events, expected after cuts, to obtain Nb ckg (cut) = 
1.8 ± 0.75 ± 0.40. From the method used in this section, 
the background is given by N bckg (cut) = e\ ckg -(l- f t t)-N 
which gives Nb ckg [cut) = 1.7 ± 0.8 ± 0.5. The two ways 
of estimating the background share in common only the 
fraction expected after A, Ht cuts. The total amount 
of background to the 4 or more jets channel is estimated 
in one case by the fit to jet multiplicities, in the other 
by a fit to the A, Ht plane. It is worth emphasizing 
that neither method uses the total number of W + 
>4 jet events predicted from Montecarlo. We can 
compute the tt cross section after doing a background 



subtraction (averaging the background obtained with 
the 2 methods) to get a t i — 7.3 ± 7.2 ± 3.3 pb for 
rritop = 160 — 180 GeV. The value should be increased 
by 25% for m top = 140 GeV. 

5. Conclusion 

Two different ways of analyzing lepton+jets final states 
lead to a similar estimate of the top production cross 
section in pp interactions at 1.8 TeV center-of-mass 
energy: 

a ti = 6.4 ± 9.8 ± 4.0 pb 

a t i = 8.1 ± 6.9 ± 3.8 pb 

and of the background after A > .05, H T > 140 GeV 
cuts: 

1.8 ±0.9, and 1.7 ±0.9 

with 4 events observed. 

The data shows a small excess of events above the 
expected background but the excess is not statistically 
significant. The excess above background corresponds 
to a top cross section of a t t — 7.3 ± 7.2 ± 3.3 pb. These 
results should be considered preliminary. 
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